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The galilean genesis scenario is an alternative to inflation in which the universe starts expanding 
from Minkowski in the asymptotic past by violating the null energy condition stably. Several 
concrete models of galilean genesis have been constructed so far within the context of galileon-type 
scalar-field theories. We give a generic, unified description of the galilean genesis scenario in terms 
of the Horndeski theory, i.e., the most general scalar-tensor theory with second-order field equations. 
In doing so we generalize the previous models to have a new parameter (denoted by a) which results 
in controlling the evolution of the Hubble rate. The background dynamics is investigated to show 
that the generalized galilean genesis solution is an attractor, similarly to the original model. We 
also study the nature of primordial perturbations in the generalized galilean genesis scenario. In all 
the models described by our generalized genesis Lagrangian, amplification of tensor perturbations 
does not occur as opposed to what happens in quasi-de Sitter inflation. We show that the spectral 
index of curvature perturbations is determined solely from the parameter a and does not depend on 
the other details of the model. In contrast to the original model, a nearly scale-invariant spectrum 
of curvature perturbations is obtained for a specific choice of a. 

PACS numbers: 98.80.Cq, 04.50.Kd 


I. INTRODUCTION 


It is fair to say that inflation [1-3] followed by a hot Big 
Bang is a standard scenario of modern cosmology. Infla¬ 
tion is attractive because the period of quasi-de Sitter 
expansion in the early universe resolves several problems 
that would otherwise indicate the need for fine-tuning. 
Moreover, curvature perturbations are naturally gener¬ 
ated from quantum fluctuations during inflation, which 
seed large-scale structure of the universe [!]. The ba¬ 
sic prediction of inflation is that the primordial curva¬ 
ture perturbations are nearly scale-invariant, adiabatic, 
and Gaussian. This is in agreement with observations of 
CMB anisotropies [5, 6]. Inflationary models also pre¬ 
dict the quantum mechanical production of gravitational 
waves [7], the detection of which would be the evidence 
for inflation. 

Despite the success of inflation, it would be reasonable 
to ask whether only inflation can be a consistent scenario 
compatible with observations. It should also be noted 
that an inflationary universe is past geodesically incom¬ 
plete [8] and so the problem of an initial singularity still 
persists. From this viewpoint, various alternative scenar¬ 
ios have been proposed so far, such as bouncing models. 
Although such models can eliminate the initial singular¬ 
ity, many of them are unfortunately plagued by insta¬ 
bilities originated from the violation of the null energy 
condition (NEC), the growth of shear, and primordial 
perturbations incompatible with observations [9]. 

In the context of cosmology, the violation of the NEC 
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implies that 


where H is the Hubble rate and t is cosmic time. This 
signals ghost instabilities in general relativity. Recently, 
however, it was noticed that in noncanonical galileon- 
type scalar-field theories the NEC can be violated sta¬ 
bly,^ and based on this idea, Creminelli et al. proposed a 
novel, stable alternative to inflation named galilean gen¬ 
esis [II]. (See also Ref. [12].) In the galilean genesis sce¬ 
nario, the universe is asymptotically Minkowski in the 
past and starts expanding from this low energy state. As 
such, this scenario is devoid of the horizon and flatness 
problems. Aspects of galilean genesis have been studied 
in Refs. [13-17] and the original model has been extended 
in Refs. [18-20] to possess improved properties. See also 
Refs. [21-29] for other interesting NEC violating cosmolo¬ 
gies in galileon-type theories and Ref. [30] for a related 
review. 

In this paper, we introduce a unified treatment of the 
galilean genesis models and give a generic Lagrangian ad¬ 
mitting the genesis solutions. This is done by using the 
Horndeski theory [31], which is the most general scalar- 
tensor theory with second-order field equations. Our gen¬ 
eralized galilean genesis Lagrangian contains four func¬ 
tional degrees of freedom and a constant parameter de¬ 
noted a. This parameter determines the behavior of the 
Hubble rate. For specific choices of those functions and 
a = 1, our Lagrangian reproduces the previous models 
explored in Refs. [11, 18-20]. As is often the case with in¬ 
flation alternatives, it turns out that the galilean genesis 


^ The NEC can be violated stability at least within linear per¬ 
turbation analysis. However, at nonlinear order, it is not clear 
whether there are no instabilities [iO]. 
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models in general fail to produce nearly scale-invariant 
curvature perturbations. We show, however, that with 
an appropriate tuning of a it is possible to have a slightly 
tilted spectrum consistent with observations. 

The Horndeski theory was developed about forty years 
ago [31] and was revived recently as the generalized 
galileon theory [32]. The equivalence of the two theories 
was proven for the first time in Ref. [33]. The action of 
the Horndeski theory is given in the generalized galileon 
form by 


S — J d'^Xyy—g (£2 + .^3 + -^4 + ^ 5 ) I (2) 

with 

C2 = G2{^,X), /:3 = -G 3 (<^,X)n<^, 

£4 = Gii^, X)R + Gax [(□<^)" - , 

£5 = G5(</>, - ^G5X [(□</>)' 

-3n<^(V^V,0)2 + 2(V^V,<^)3] , (3) 

where R is the Ricci scalar, is the Einstein tensor, 
and each Gi(^, X) {i = 2, 3,4, 5) is an arbitrary function 
of the scalar field ^ and X := —g^''dfj_(l)d,^(j)/2. We use 
the notation Gix to denote dGijdX. 

The plan of this paper is as follows. In the next section, 
we present a generic Lagrangian that admits the gener¬ 
alized galilean genesis solution. In Sec. Ill we analyze 
the background evolution analytically and numerically 
and show that the generalized galilean genesis solution is 
the dynamical attractor for a wide range of initial con¬ 
ditions. Primordial tensor and scalar perturbations from 
generalized galilean genesis are studied in Sec. IV, and 
the curvaton mechanism in the genesis scenario is briefly 
discussed in Sec. V. We draw our conclusions in Sec. VI. 


II. GENERALIZED GENESIS SOLUTIONS 


In Ref. [34] it was noticed that the genesis solution (5) 
is obtained generically in the subclass of the Horndeski 
theory with 

G2=e^^^g2iY), G3 = e^^^gsiY), 

G4 = ^+e2^^54(y), G5=e-^^^g,{Y), (6) 

where each gi {i = 2, 3,4, 5) is an arbitrary function of 
Y := (7) 


This extends the Lagrangian given in Ref. [16] to include 
the Horndeski functions G 4 and G 5 . The Lagrangian (4) 
and the DBI conformal galileon theory are included in 
the general framework defined by ( 6 ) as specific cases. 

In this paper, we further generalize ( 6 ) and consider 

G 2 = G 3 = e^‘^^^g3iY), 

G 4 = M + e^‘^^^g4Y), Gs = e-^^*g^{Y), ( 8 ) 


where a (> 0) is a new dimensionless parameter. The 
four functions, 32 , ffa, 54 , and 55, are arbitrary as long as 
several conditions presented in this section and in Sec. IV 
are satisfied. We assume, however, that 54(0) = 0, so that 
G 4 —Mpj/2 as y —0. The Horndeski theory with ( 8 ) 
admits the following generalized galilean genesis solution: 


gA 0 


1 1 TT ^ ^0 


(—00 <t < 0 ),( 9 ) 


for large [t], where Yq and Iiq are positive constants. We 
see that V ~ Ig for this background. The parameter a 
in the Lagrangian results in controlling the evolution of 
the Hubble rate. The scale factor is given by 


a — 1 -j- 


1 hp 
2a 


( 10 ) 


The original model of galilean genesis is constructed 
by using the Lagrangian of the form [11, 18] 

]Vf 2 

£ = + fie^^'^X + f2X^ + faXDc^, (4) 


where /i, / 2 , fs, and A are constants. (We have changed 
notations of Refs. [11, 18].) The above Lagrangian has 
the genesis solution. 




const 

^T’ 




const 


(—00 < t < 0), 


( 5 ) 


for large jtj. (We have a degree of freedom to shift the 
origin of time: t —>■ t — to.) The scale factor is given 
by a ~ 1 - 1 - const/(—t)^, describing the universe that 
starts expanding from singularity-free Minkowski in the 
asymptotic past. The same genesis solution can also be 
obtained from the DBI conformal galileons [19, 20]. 


and hence the solution describes the universe that starts 
expanding from Minkowski in the asymptotic past, sim¬ 
ilarly to the original galilean genesis solution which cor¬ 
responds to the case of a = 1. The “slow-expansion” 
model considered in Ref. [35] is reproduced by taking the 
particular functions gi with a = 2. We thus obtain a one- 
parameter family of the generalized genesis solutions as 
an alternative to inflation. Note that, although the evo¬ 
lution of the scale factor is very different from quasi-de 
Sitter, the universe in this scenario is also accelerating: 
dt{aH) > 0, and hence fluctuation modes will leave the 
horizon during the genesis phase. 

Substituting Eq. (9) to the background equa¬ 
tions (A1)-(A3) and picking up the dominant terms at 
large [t], we have 

£ ~ e 2 (“+i)^‘^p(yo) 0 , ( 11 ) 

V-2giYo)H + e^^‘^+^^^^piYo)~0, ( 12 ) 
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where 


III. BACKGROUND EVOLUTION 


P{Y) 

:= 2Yg'^ - 32 - 4Ar {ag^ - Yg'^) , 

(13) 

P{Y) 

:= 32 - 4 aAT 33 



+S{2a+l)\^Y{agi-Yg'^), 

(14) 

G{Y) 

:= M|i-4Ay(35 + T 3 '), 

(15) 


To see whether or not the generalized genesis solu¬ 
tion presented in the previous section is an attractor, 
we trace the background evolution starting from generic 
initial conditions. 


an overdot stands for differentiation with respect to t, 
and a prime for differentiation with respect to Y. The 
constant Yq is determined as a root of 


piYo) = 0 , 

and then Hq is determined from Eq. (12) as 

. ^_ 1 PjYo) 

° 2{2a + l){2X^Yo)^+^g{Yo)' 


(16) 


(17) 


As will be seen shortly, this background is stable for 
G{Yo) > 0. Therefore, the above NEC violating solution 
is possible provided that 


A. Analytic argument 

Let us begin with a simplified discussion neglecting 
gravity, i.e., the effect of the cosmic expansion [16]. It is 
convenient to introduce a new variable 

^ := e-^'^ (> 0 ). ( 21 ) 

In terms of ip we have Y = pp'/{2}?). Eor any homoge¬ 
neous solutions the scalar-field equation of motion (A4) 
with the functions ( 8 ) can be written as 


p{Yo) < 0 . (18) 

As will be demonstrated in the next section, the gener¬ 
alized genesis solution will develop a singularity ^ oo 
at some t = tsing, as in the original genesis model. 
We therefore assume that the genesis phase is matched 
onto the standard radiation-dominated universe before 
t = tsing, ignoring for the moment the detail of the reheat¬ 
ing process. In conventional general relativity, matching 
two different phases can be done by imposing that the 
Hubble parameter is continuous across the two phases. 
However, the matching conditions are modified in general 
scalar-tensor theories as second-derivatives of the metric 
and the scalar field are mixed in the field equations. The 
modified matching condition [34] reads 

J2a+l)\<l> rYo 

= g{Yo)H -^ ^/^g3{y)dy 

+2A0e2“^<^ {ag4 - Toff^), (19) 

and we require that the subsequent radiation-dominated 
universe is expanding: iLrad > 0. This condition trans¬ 
lates to 

r^o 

- 32 - 2 Ayo <?3 + ( 2 a + l)Av^ / -^dy > 0 . ( 20 ) 

do \/v 

It is easy to see that in the case of a = 1 all the expres¬ 
sions presented above reproduce the previous results [34] . 

Before closing this section, let us emphasize that (gen¬ 
eralized) galilean genesis has the Minkowski phase only 
in the asymptotic past. The true Minkowski spacetime 
solution corresponds to the special case of T = 0 , i.e., 
(j) = const. The Y = 0 solution is found only if 32 ( 0 ) = 0. 
One may wonder if the true Minkowski vacuum {Y = 0) 
in our neighborhood begins to expand to form a genesis 
universe (T = Iq > 0 )- This is forbidden because the 
two different stable solutions cannot be interpolated, as 
argued in Ref. [16]. (See, however, Ref. [17].) 


— U- 2 (“ 


dt . 


iIj-P°‘+P p{Y) 


= 0 . 


Integrating this, we obtain 
p{Y) = 


( 22 ) 


(23) 


where C is an integration constant. Equation (23) de¬ 
fines a curve in the {ip, ip) space for each C, as shown in 
Fig. 1. With an initial condition {ipi,ipi) away from the 
genesis solution, the integration constant is determined 
as C = ip^ ^^“~''^^/3('!/i?/ 2A^). If 'i/’ < 0 initially, the scalar 
field rolls along the curve toward ip —>■ 0 , i.e., p —>■ 0 . 
Hence, this solution approach to one of the genesis solu¬ 
tions which are denoted as horizontal lines {ip = const) 
in the {ip, ip) plane. If ■i/' > 0 initially, the scalar field rolls 
the opposite way along the curve and goes further away 
from the genesis solutions. This is the time reversal of 
the Ip < 0 solutions. 

The above analytic argument implies that the genesis 
solution is the attractor for initial conditions such that 
Ip < 0 (<J4> {e^^y > 0). In the next subsection we perform 
numerical calculations to show that this is basically true 
even if one takes into account of the effect of gravity. The 
numerical analysis also allows us to see the final fate of 
the genesis solutions for which the effect of the cosmic 
expansion cannot be ignored. 


B. Full numerical analysis 

In the Horndeski theory with ( 8 ) the Eriedmann equa¬ 
tion can be written as 

£ = e2(“+i)^‘^/5(y)-H6i7(^e2“^‘^ci(r) 

-3ij2 [c 2{Y) + e^°‘^’^d2{Y)] + 2H^Pe-^^*C3{Y) 
= 0, (24) 






4 


ij) 



FIG. 1: Examples of the curves defined by Eq. (23). Hori¬ 
zontal dashed lines correspond to the genesis solutions. 

where 

Cl = Yg'^-2aXg4 + 2{3-2a)XYg'^ + 4:XY^g'l,(25) 

C 2 = M^,-l2XYg5-28XY^g'^-8XY^g'', (26) 

C 3 = 5Yg', + 2Y^g'', (27) 

d2 = 2g4-8Yg'^-8Y^g''. (28) 

Equation (24) is exact and hence can be used even if 
the background evolution is away from the genesis solu¬ 
tion. Similarly, one can substitute Eq. ( 8 ) to the evolu¬ 
tion equation V = 0 and the scalar-field equation of mo¬ 
tion to write straightforwardly the exact equations for the 
background. The resultant equations are integrated nu¬ 
merically, giving the background evolution starting from 
generic initial conditions. 

Given the initial conditions {(j}{tQ), (j)(to)), the initial 
value for H is determined from the Friedmann equa¬ 
tion (24). Therefore, the initial values {4>{to), 0(to)) must 
be chosen in such a way that Eq. (24) admits a real root 
H. Equation (24) is quadratic in if ^5 = 0 and cubic 
if 95 7 ^ 0 - In both cases, the discriminant V for 1 

is given by 

V = e2(“+i)^’^c2(E)/3(y) -h C)(e2(2“+i)Ad)^ ( 29 ) 

In the 35 = 0 case, the initial data {(j){to), (p{tQ)) must lie 
in the region where 2? > 0 is satisfied. In the gs 0 case, 
the Friedmann equation has at least one real root for any 

Concrete numerical examples are presented in Figs. 2- 
4. In Figs. 2 and 3 we show the cases where the Fried¬ 
mann equation is quadratic in H. The shaded regions 
{T> < 0) cannot be accessed because H would be imagi¬ 
nary there. In Fig. 2 we have one genesis solution, while 
we have two in Fig. 3. In both cases, generalized galilean 





FIG. 2: Numerical results of the background evolution for the 
model with 32 = —Y + Y'^, gs = Y, and qa = gs = 0. The 
parameters are given by Mpi = 1, A = 1, and a = 1. 



* 

FIG. 3: Numerical results of the background evolution for the 
model gi = —Y 3Y^ — Y^, g^ = Y, and g 4 = gs = 0. The 
parameters are given by Mpi = 1, A = 1, and a = 2. 

genesis is the attractor for ■0 < 0. At late times where 
0 1 , the numerical solutions are no longer approxi¬ 

mated by Eq. (9), and within a finite time the Hubble 
rate H diverges. As the Friedmann equation is quadratic 
in H, we have two branches of the solutions, one of which 
may be contracting initially {H < 0). The cosmological 
evolution nevertheless approaches the same genesis solu¬ 
tion and the trajectories in the ( 0 , tp) space are almost 
indistinguishable. 

The behavior of the models with 55 0 0 is more com¬ 
plicated, as illustrated in Fig. 4. In the white region, we 
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FIG. 4: Numerical results of the background evolution for 
g 2 = —Y +y^, gs = Y,g 4 = 0, and g^ = —Y. The parameters 
are given by Mpi = 1, A = 1, and a = 1. 


have T) > 0 and so there are three possible choices for the 
initial value of H. Two of the three branches converge 
to the genesis solution similarly to the 35 = 0 case, as 
shown as the black lines in Fig. 4. Also in this case we 
find H ^ 00 within a finite time. However, the remain¬ 
ing one branch never converges to the genesis solution. 
The corresponding examples are shown as the gray lines 
in Fig. 4. In the shaded region, we have V < 0 and there 
is only one possible initial value for H at each point, 
which corresponds to the latter branch. Therefore, the 
generalized galilean genesis solution can be a dynamical 
attractor for the initial data in the white {V > 0) region. 

We thus conclude that the galilean genesis solution is 
the attractor provided that ^ < 0 (<f4> > 0 ) ini¬ 

tially, though the situation in the presence of is in¬ 
volved. In the inflationary scenario, usually it does not 
matter which direction the scalar field rolls initially, but 
the universe must be expanding initially. In contrast to 
the case of inflation, the galilean genesis scenario allows 
both for expanding and contracting universes at the ini¬ 
tial moment, while the scalar field must roll in a particu¬ 
lar direction initially. As far as we have investigated nu¬ 
merically, all the solutions develop a singularity H ^ 00 
at some time t = tging in the future. In passing, we have 
checked that the numerical examples in Figs. 2 and 4 sat¬ 
isfy the stability conditions presented in the next section. 


C. Spatial curvature 

We have so far neglected the spatial curvature. In this 
subsection, let us justify this assumption by showing that 
the spatial curvature does not interfere with the evolu¬ 
tion of the genesis background. We will use the cosmo¬ 


logical background equations with the spatial curvature 
K in the Horndeski theory, which are summarized in the 
Appendix A. 

Let us take an initial condition such that H is suffi¬ 
ciently small in the equation of motion for (j) and ' > 
0. Then, in a universe with K ^ 0 the equation of motion 
for (j) can be written as 


- ^‘2a\4> ^4{Y) _ K.^{Y) 

dt a? 


= 0, (30) 


where 


lCi{Y) := 6{g4-2Yg'4)K, (31) 

/C 5 (T) := -12XY{g5 + Yg',)K. (32) 


Even if p ^ /C 5 at the initial moment, 

the curvature terms become smaller relative to the p term 
as the scalar field rolls. Thus, we have the same attractor 
solution F = Fq satisfying piY^) = 0 , i.e., ~ 

Along this attractor, the evolution equation reads 


2giYo)H + e2(“+i)^^p(Fo) + + iXYog,{Yo)] J ~ 0, 

(33) 


where we assumed that H ^ H^. Equation (33) im¬ 
plies that the curvature term becomes subdominant as 
the scalar field rolls, and as a result H is determined by 
the p term, recovering the evolution of the genesis back¬ 
ground. Thus, the flatness problem is resolved in the 
genesis model. 


D. Anisotropy 


In conventional cosmology, an initial anisotropy is 
wiped out during inflation [36]. However, in alternative 
scenarios such as bouncing cosmology, it is often prob¬ 
lematic that the initial anisotropy grows in a contract¬ 
ing phase [37] (see however [26]). In this subsection we 
will show that adding the initial anisotropy on the gen¬ 
eralized galilean genesis solution does not destabilize the 
background evolution. 

We consider the Kasner metric 


ds^ = —dt^ 


A«id)da;2 + e2^^(*)dy 


2 , .2fl3(t)d^2 


(34) 


where it is convenient to write 

01 =/3+ + y3/3_, 02 = /3+-^/3-, 03 =-2/3+.(35) 


If the deviations from the genesis background are not 
large, it follows from Eqs. (B4) and (B5) that 


'gp+ - 2e-^^*^Yog’, (^1 - pl) 

^ [a/3-+4e-2^^(^Fog'/3+/3_ 


0, (36) 

0. (37) 
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In the models with g'^ = 0, this simply gives 


B. Scalar perturbations 


/3+, /j_ ^ const, (38) 

so that the initial anisotropy dilutes as 9i ~ (—t)- In the 
models with g'^ ^ 0, we have the following possibilities: 

(/3+,/3_) = (0,0), (&, 0), i-^b, ^6), {-h, -^6),(39) 
where 


b 


g 




(40) 


In this case, for nonzero /3± the initial anisotropy can 
grow logarithmically: 6i ~ ln(—t). However, this should 
be compared with Ina ^ (—1)“^“; we see that the log¬ 
arithmic growth of 9i does not spoil the genesis back¬ 
ground. 


IV. PRIMORIDAL PERTURBATIONS 

Let us now study the behavior of primordial tensor 
and scalar perturbations around the generalized genesis 
background to obtain predictions of our scenario as well 
as to impose stability conditions. To do so, we utilize the 
general quadratic action for cosmological perturbations 
in the Horndeski theory derived in Ref. [33]. 


A. Tensor perturbations 


The quadratic action for tensor perturbations hij in 
the genesis phase is given by 



dt9?xa^g{Yo) 




(41) 


where 


2 _ Mgi + 4Ayog5(yo) . 

‘ g{Yo) ^ ^ 

and note that a ~ 1. It can be seen that stability against 
tensor perturbations is assured if 


g{Yo) > 0, (43) 

M2i + 4AFoff5(Lb) >0, (44) 


are satisfied. 

Since both t/(Lb) and are constant during the gene¬ 
sis phase, the tensor perturbations are effectively living in 
Minkowski without regard to a and the concrete form of 
gi{Y), and consequently amplification of quantum fluctu¬ 
ations does not occur as opposed to the case of quasi-de 
Sitter inflation. This means that no detectable primor¬ 
dial gravitational waves are generated from our generic 
class of the genesis models. 


The quadratic action for the curvature perturbation ( 
in the unitary gauge is given by 



dtd^xa^gs 


c" 



(45) 


where with some manipulation C/s and in the genesis 
phase are written as 


gs 


with 


(2a + l)Xe{Yo) 


YoCiYo) 

C(Yo)p(Yo) 

aYo)p'{Yoy 


p'(ro)e- 2 “^^ 


(46) 

(47) 


:= 


Yg{Y) 

p{Y) 


(48) 


Equations (46) and (47) show that oc (—1)^“ and = 
const. It follows from Eqs. (18) and (43) that ^(Lb) > 0. 
We thus find that stability against scalar perturbations 
is guaranteed if 


p'iYo) > 0, (49) 

e{Yo) < 0, (50) 


are fulfilled. We can choose the functional degrees of 
freedom so that this is possible. 

Let us evaluate the power spectrum of C. To simplify 
the notation, it is convenient to write t/g = Al(—<)^“, 
where .4 is a constant deduced from Eq. (46), the value 
of which depends on the model, i.e., a and the concrete 
form of gi(Y). The equation of motion derived from the 
action (45) is given by 

2 a • 

Cfc H—= 0, (51) 

where we moved to the Fourier space. This equation can 
be solved to give 


= I - a. (52) 

where is the Hankel function of the first kind and 
the positive frequency modes have been chosen. On large 
scales, jcsfctj 1, we have 


Cfc — Alfc + Bk{—t) 


l-2a 


(53) 


where 
Ak : 
Bk : 


= -^ 2 " 





2A TT 


{csk)-y 


(54) 


_ n — — 1 I 

2A 

x{csky 


i cos(7ri2)r(—I/) 


r(i. + i) 


(55) 
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If 0 < Of < 1/2, the second term in Eq. (53) decays as is 
common to usual cosmologies, leaving the constant mode 
at late times. Thus, in this case the power spectrum is 
given by 


Vcik) 


22i/-4^-2i/ 




3-2i^ 




(56) 


and the spectral index is found to be 


Us = 2a + 3, 


(57) 


yielding a blue spectrum incompatible with observations. 

The case of a > 1/2 is more subtle, because the second 
term in Eq. (53) grows and dominates on large scales. 
This is what happens in the original galilean genesis 
model (a = 1) [11]. To extract the late-time amplitude of 
Cfe, let us consider the following situation. Suppose that 
the genesis phase terminates at t = tend and is matched 
onto some other phase. We assume that the scalar field 
is homogeneous on the t = tend hypersurface. In the sub¬ 
sequent phase, the curvature perturbation on large scales 
may be written as 

/■°° dt' 

a = Ck-Dkj^ ^^(tOo^M, 


where we do not specify Gsit) for t > tend, but assume 
that Eq. (58) gives the constant and decaying modes 
and hence the integral converges. The late-time ampli¬ 
tude is given by Ck ■ The matching conditions [34] imply 
that and GsCk are continuous across the two phases 
(cf. [38]). It is then straightforward to obtain Ck = 
Afe + Bfc(-tend)^-"“(l+T) ^ Bfe(-te„d)l-2“(l+T), 
where 


I := {2a 



a^(tend)-4(-tend)^° dt' 
a?‘{t')Gs{t') [tend] 


(59) 


is independent of k. We may thus use the estimate 


Vdk)^CxVc{k)\,^,^^^, (60) 

with C being some ^-independent factor. The power spec¬ 
trum evaluated at the end of the genesis phase is given 

by 




2-2'^-4cf r 2 (-Z/) 
■n'^A 


^end I 


\iv 1,3+21/ 


(61) 


so that 


Us = 5 — 2a. (62) 

Although the overall amplitude depends on the details 
of the model construction, the spectral index depends 
only on a and not on the concrete form of gi{Y). We 
have an exactly scale-invariant spectrum for a = 2, and 
this is in sharp contrast to the original galilean genesis 
model having a = 1, which produces a blue-tilted spec¬ 
trum of curvature perturbations. A particular realization 
of a = 2 is found in Ref. [35, 39], where the same con¬ 
clusion is reached. Taking a = 2.02, one can obtain the 
nearly scale-invariant, but slightly red-tilted, spectrum 
with Us ~ 0.96. 


V. CURVATON 


In the previous section we have seen that the nearly 
scale-invariant spectrum for curvature perturbations is 
possible only in the case of a ~ 2. In the other cases we 
need to consider an alternative mechanism such as the 
curvaton in order to obtain a scale-invariant spectrum. In 
this section, we study slightly in more detail the curvaton 
coupled to a conformal metric, the basic idea of which 
was proposed earlier in Ref. [11]. A similar mechanism 
was proposed in Ref. [14]. 

To make a scale-invariant power spectrum, we intro¬ 
duce a curvaton field a coupled to the conformal metric, 

(63) 


where /3 is a constant parameter which is assumed to be 
close to unity, ~ 1. Assuming the simplest potential, 
we consider the following action for cr: 


Act = / 


-^g^^d^ad^a - 


• (64) 


The conformal metric (63) implies that the effective scale 
factor for the curvaton is ^ {—t)~^ with /3 ~ 1, so 
that a lives effectively in a quasi-de Sitter spacetime. 

The equations of motion for the homogeneous part a = 
CTo(t) is given by 


do + (2/3A(/)-|-3i7)(To + e^^''''^m2crQ = 0. (65) 

On the genesis background, one can ignore H ~ 
(—t)“(2a+i) relative to At/ ~ {—t)~^, leading to 

2/3 , m? , . 

“ T"" + P+lKHjp-"" = "■ <“> 

The effective Hubble rate for the curvaton is ~ \y^2Yo. 
Eor the “light” curvaton with 

< A^Fo, (67) 


we thus have ctq — const and the other independent so¬ 
lution decays quickly. 

The energy density and pressure of a are given by 

P. = + ( 68 ) 

P. = (69) 

Equations (11) and (12) imply that the dominant part 
of the cosmological background equations grows as ^ 
(—t)“2(“+i). Thus, in order for the (initially subdom¬ 
inant) curvaton not to spoil the genesis background as 
time proceeds, we require that 


a + l>2l3. (70) 

The fluctuation of the curvaton, Sa{t,x.), obeys 


Sa - —Sa — V^Jcr 

t 


[Ay2l/(-t)]2/3 


Sa = 0. (71) 










Neglecting the mass term, this can be solved in the 
Fourier space to give 

Sau = ^ (72) 


where the positive frequency modes have been chosen. 
Thus, the power spectrum of the curvaton fluctuations is 


Vsa{k) 


7r3 


(73) 


and we find 


n, = i- 2/3. (74) 

In the case of /3 = 1, the effective scale factor for the 
curvaton is that of exact de Sitter, and hence the power 
spectrum is exactly scale-invariant, as is expected. Tak¬ 
ing (3 = 1.04 we obtain Ug = 0.96. The curvaton 
fluctuations can be converted into adiabatic ones after 
the genesis phase, where a behaves as a conventional 
scalar held in a true expanding universe, in the same 
way as the usual curvaton held in the inhationary sce¬ 
narios. Note, however, that due to the restriction (70) 
the present curvaton mechanism works only for the mod¬ 
els with a > 2 — Us {> 1). 


VI. CONCLUSIONS 

In this paper, we have extended the galilean genesis 
models [11, 18-20] and constructed a generic Lagrangian 
from the Horndeski theory that admits the generalized 
galilean genesis solution. In generalized galilean gene¬ 
sis, the universe starts expanding from Minkowski in a 
singularity free manner with the increasing Hubble rate 
H ^ (—t)“^2“-i-i)^ where a (> 0) is a new constant pa¬ 
rameter in the Lagrangian. We have investigated the 
background evolution and shown that the generalized 
galilean genesis solution is the attractor for a wide range 


of initial conditions. In particular, we have seen that the 
spatial curvature and an initial anisotropy do not hinder 
the evolution of the genesis phase. 

We have then studied the primordial perturbations 
from the generalized galilean genesis models. From the 
quadratic actions for cosmological perturbations we have 
imposed several stability conditions on the functions in 
our generic Lagrangian. In contrast to the case of quasi- 
de Sitter inflation, tensor fluctuations are not amplified 
in the genesis phase in all the galilean genesis models 
we have constructed, and hence no detectable primor¬ 
dial gravitational waves are expected. The evolution of 
the curvature perturbation (/ depends on the parameter 
a and has turned out to be more interesting. In the 
case of a > 1/2, C grows on large scales, as in the orig¬ 
inal galilean genesis model {a = 1) [11]. The tilt of the 
power spectrum at the end of the genesis phase is given 
by = 5 — 2a, irrespective of the other details of the 
model. Thus, we have a slightly red-tilted spectrum for 
a >2. In the case of a < 1/2, the constant mode domi¬ 
nates on large scales as in conventional cosmology. In this 
case, the power spectrum has been shown to be always 
blue-tilted. We have also discussed the possibility of the 
curvaton mechanism in the generalized galilean genesis 
scenario. 

We have ignored the reheating process in our scenario. 
It would be interesting to explore how the universe re¬ 
heats and how matter is created at the end of generalized 
galilean genesis. Since the Lagrangian defined by (8) ex¬ 
cludes a cosmological constant, it is not clear how the 
genesis phase is connected finally to the late-time uni¬ 
verse described by the ACDM model. These are the open 
questions. 
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I 

Appendix A: Cosmological background equations in the Horndeski theory 


The cosmological background equations in the Horndeski theory are given by [33] 

f = 0, V = 0, 


(Al) 


where 

£ := 2XG2X -G 2 + dX^HGsx - 2XG3^ - QH^Gi + 2AH‘^X{Gix + XGixx) - l2HXij)Gi4,x - QH^Gi^ 

+2H^X^ {bG^x + 2XG5 xx) - (SGg^ + 2XG5^x), (A2) 

r := G 2 - 2 X (Gs^ + ^Gsx) + 2 (sH^ + 2 H^ G 4 - l2H^XGix - '3HXGix - 8 HXG 4 X - SHXXG^xx 

+2 -f 2H(/>'j G40 -b 4AG400 -b 4Y ((/ - 2iJ0) G4^x - 2X {2H^^ -b 2HH^ + G^x 

-AH^X^G^xx + 4HA (a - HX^ G^^x + 2 [2 {HXy+ 3H^X] Gg^ + diJA^Gg^^. (A3) 




9 


The equation of motion for (j) takes the form 


1 d 
dt 


(a^j) = 


where 


and 


J := <j)G 2 x + 6 HXG 3 X - 2^G3^ + {G 4 X + 2XG4xx) - UHXGa^x 
+2H^X (3G5X + 2XG5xx) - {G^^ + XG^^x ), 

P(j, '■= G 20 — 2X ^G300 + (fiG^fjjx^ + 6 { 211 ^ + Gifj, 

+6H (^X + 2HX^ GicfiX — QH^XG^^^ + 2H^ X^G^^x- 


(A4) 


(A5) 


(A6) 


The above equations are for the spatially flat background. In the open {K = —1) and closed {K = 1) cases, the 
corresponding equations are given by 


and 


where 


J__d 

^dt 


= f + A£ = „, 


6K • 

’3 j+ {cj,G4x + HXG5x-^G5^) 


— Pa + 


SKdPx 

a? d(j) ’ 


(A7) 

(AS) 


Gt := 2 [g4 - 2AG4jf - A (p^Ggx - G 50 )] , (A9) 

Pt ■■= 2 [g4-a(^G5;c+G 5^)] . (AlO) 

The cosmological background equations for the open and closed models are derived for the first time in this paper. 


Appendix B: Anisotropic Kasner universe in the Horndeski theory 


We derive the basic equations governing the evolution of an anisotropic Kasner universe in the Horndeksi theory. 
We consider the following metric: 

ds2 = ^ ^2 jg2(/3+ + V3/3_)^^2 ^ g2(/3+- V3/3_ )^y2 ^ g-4/3+d^2 

Substituting this to the Horndeski action, we obtain 

S — S'lso T ‘S'aniso j 

where S'lso is identical to the action for the homogeneous and isotropic metric and Saniso is given by 

N 


‘5'aniso — / X 


G 4 - 2AG4X - ^AG5x + XG^a) (/3+ + P-) - ^X^G^x - S/S+Z?" ) 


Note that X here should be understood as A = <jp/2N‘^. 

Varying the above action with respect to j3± and setting A = 1, we obtain 

^ {a^ [gTP+ - 2X^G3X (pl-Pi)]] = 0, 


^ \g^p_ + AXpG^xP+P-] } = 0. 


(Bl) 

(B2) 

(B3) 


(B4) 

(B5) 


Varying the action (B2) with respect to A, a, and (j), one can also derive the cosmological evolution equations with 
shear contributions. I_ 
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